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Abstract 



' We report on QCD radiative corrections to the processes pp — > W + 2 jets and — > Z + 2 jets 

C ■ at the Tevatron. These processes are included in the Monte Carlo program MCFM, which allows 

the calculation of any infra-red finite variable at next-to- leading order. Due to a better theoretical 



description of jets at next-to-leading order, some distributions exhibit significant corrections. As 



^ ' expected, the unphysical dependence of theoretical predictions upon the renormalization and fac- 

a: 

torization scales is greatly reduced compared to leading order. As an example of the predictions 
that may now be made with MCFM, we present a next-to-leading order estimate of the heavy 
flavor content of jets produced in association with vector bosons. 

H 

03 1 PACS numbers: 13.38 -b, 12.38 -t 



Electronic address: johnmc@hep.anl.gov 



^Electronic address: gllis@fnal.g- 



I. INTRODUCTION 



In this paper we report on the results of a calculation of the next-to-leading order QCD 
corrections to the processes, 

p + p — > W + 2 jets, 

p + p -> Z/^ + 2 jets. (1) 

These reactions will be investigated at the Tevatron, i.e. pp collisions at \/s = 2 TeV. 
Our calculations are equally applicable to the LHC (pp collisions at y/s = 14 TeV). We 
plan to consider the LHC in a subsequent publication. The results are obtained from the 
Monte Carlo program MCFM which allows us to obtain full predictions for any infra-red 
safe variable. In order to obtain fully differential distributions, to which experimental cuts 
may be applied, various decay modes of the Z/j* intermediate states are included, 

Z/ 7 * -> e~e + 
Z/Y -> bb 

z -> y~]*w (2) 

i 

as well as, 



W + -> z/ e e+ (3) 

for the W + decay. In this paper we will only report on leptonic decays of the vector bosons. 
We use the approximation of massless leptons, so that our results are also valid for the 
decays W + — > W + — > u t t + in this approach. 

Because of their phenomenological importance, processes involving the production of 
vector bosons and jets have been considered by many authors. W-hoson production with 
two jets was considered at leading order in ref. |l|, @, [3|. The same process involving jets 
at large rapidity was considered in ref. |||. Vector boson production in association with 
n-jets for n < 4 is calculated at leading order in refs. || ||. In refs. f7|, j|, |^, H^] predictions 



were made for processes involving a vector boson recoiling against one jet at next-to-leading 
order. In refs. [11, [12[ predictions were made for processes involving W bosons and one 



heavy quark at next-to-leading order. However to the best of our knowledge this is the first 
paper to calculate vector boson processes with two jets at next-to-leading order. 



In performing these calculations we have used the results of other authors for the crossed 
reactions e + e~ — > 4 partons and e + e~ — > 5 partons ||14|| . Even with the amplitudes in 



hand, the implementation in a Monte Carlo program requires considerable effort. 

In order to highlight the similarities between the effects of radiative corrections on the 
W/Z + 2 jet rate and the W/Z + 1 jet rate, we will also present some results for the 
latter process. Such corrections have been known for some time |Tt| and have provided an 



invaluable tool for studies at the Tevatron. The inclusion of the W/Z + 1 jet processes in 
MCFM was useful to understand the issues to be faced in implementing the more complicated 
W/Z + 2 jet processes. 



We will also tie together our results with previous predictions made using MCFM [|T5| , |16 
to provide a consistent next-to-leading order prediction for the heavy flavor content of jets 
produced in association with a W/Z. Specifically we consider vector boson events containing 
two tagged 6-jets. This quantity will be used in assessing the backgrounds to a number of 
new physics searches at the Tevatron [17| . Experimental studies [18] have so far relied upon 



leading order predictions as theoretical input. 

II. Hadron + Hadron -> W/Z + 2 jets 

A representative sample of the Born diagrams for the processes, 

parton + parton — > W/Z/j* + 2 partons, (4) 

is shown in Fig. [TJ. In the coding of these processes into MCFM we have made an artificial 
separation between processes involving two quarks, Fig. p] (a-c) and processes involving four 
quarks Fig. [I] (d-f). This separation is motivated by the relative size of the contributions at 
leading order, illustrated in Fig. |], as well as by the relative complexity of evaluating the 
different matrix elements. The 2-quark process at leading N dominates the total, a trend 
that we find is preserved at next-to-leading order. 

We have not recalculated the virtual corrections to the basic Born process which are given 



by Bern, Dixon and Kosower in ref . |T3| . These amplitudes have been calculated in the four 
dimensional helicity scheme, which we consistently use throughout this program. 
The real corrections to the basic Born processes, i.e. the processes, 

parton + parton — > W/Z/'y* + 3 partons (5) 
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FIG. 1: Sample diagrams for the process parton + parton — * W/Z + 2 partons. As usual the vector 
boson is denoted by a wavy line. 
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FIG. 2: Color structure of the W+2 jet cross section vs. the dijet mass Mjj at leading order. 
Proceeding from the top the five curves are, the total LO result, the two quark process at O(N), 
the four quark process at O(l), the four quark process at 0(1/N) and the two quark process at 
0(1/N), where N is the number of colors. 
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FIG. 3: Sample diagrams for the process parton + parton — > W/Z + 3 partons. 



have been published in fll4| , [19| , [20| and a representative sample of the contributing diagrams 
is shown in Fig. 0. These matrix elements are incorporated in MCFM using the subtraction 



method of Ellis, Ross and Terrano |21], [22| . In this method, one constructs counterterms 
having the same singularity structure as the real emission matrix elements, starting from 
the eikonal formula for soft emission. The basic operating procedure is to create subtraction 
terms which contain the same singularity structure as the lowest order diagrams, but are 
simple enough that they are can be integrated over the phase space of the unobserved parton. 
In order to improve the cancellation, the kinematics of the counterterms are included using 
the prescription of Catani and Seymour [p3 ]. For the case of W + 2 jets there are 24 different 
counterterms, each with its own kinematic structure, (a counter event) . The kinematics of 
the final state partons in the counterevents coincide with the kinematics of the event in the 
appropriate soft and collinear limits. 
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A. Numerical checks 



The matrix elements for hadroproduction of W/Z+3 partons contain many singularities 
which are subtracted by (dipole) counterterms. Although the enumeration of these coun- 
terterms is in principle straightforward, the success of the whole program depends on it 
being implemented correctly. We shall therefore present a few details of the checks which 
we performed. 

The real matrix elements were taken from ref. \TA\ supplemented in certain cases by our 
own calculations. The coding of these matrix elements was checked by comparison with 
routines generated by the program MADGRAPH powered by HELAS []25|| . We cannot 
use the routines generated by MADGRAPH directly because the resulting code is too slow 
to implement in a Monte Carlo program which requires many calls to the matrix element 
routine. 

The next step is to verify that the numerical value of the counterterms is in fact equal in 
magnitude to the real matrix element in the singular limit. This is done by generating sets of 
points which lie in all of the potentially singular regions and then checking the cancellation 
of the event and the appropriate counterevent in each limit. 

Finally, we must add back the counterterms, suitably integrated over the phase space of 
the emitted parton. Here it is clearly important to add back exactly what has been sub- 
tracted in the previous step. We have tried to structure the code so that the comparison 
between the two steps is transparent. The integral of the counterterms over the emitted par- 
ton sub-space contains singularities, which are regulated using dimensional regularization, 
in addition to finite contributions. The code is structured so that the finite parts are closely 
associated with the singular parts. Thus the cancellation of the singular parts, (with the 
sum of the singular parts of the virtual matrix elements and the Altarelli-Parisi factorization 
counterterms for the parton distributions), provides some assurance that the finite terms are 
included correctly. 
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Parameter 


Default value 


Parameter 


Default value 


M z 


91.187 GeV 


a(M z ) 


1/128.89 


T z 


2.49 GeV 




1.16639xl0~ 5 


M w 


80.41 GeV 


gl 


0.42662 (calculated) 


T w 


2.06 GeV 


sin 2 6 W 


0.23012 (calculated) 



TABLE I: Default parameters in the program MCFM. 
III. MONTE CARLO RESULTS 



A. Input parameters 



MCFM has a number of default electroweak parameters which we use throughout this 
paper. They are given in Table |. As noted in the table, some parameters are calculated 
using the effective field theory approach [f2(f| . 



G 

e 2 = Aira(Mz), gl = 8M^^, sin0,„ = ^. 



72' 



I w 



(6) 



For simplicity we have taken the CKM matrix to be diagonal in the W + 2 jets process. As 
a consequence there are, for example, no us initial states for this case. This approximation 
is not expected to influence any anticipated analyses. For the other processes we retain only 
the Cabibbo sector of the CKM matrix, 



V ( 



CKM 



\ 



V 



(7) 



0.975 0.222 
0.222 0.975 
II 

The value of as(Mz) is not adjustable; it is determined by the chosen parton distribution. 
A collection of modern parton distribution functions is included with MCFM, but here we 



concentrate only on one of the MRST2001 [27] sets with as(Mz) = 0.119. We refer to this 
set as MRS0119, which is the label used in our program. 
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B. Basic cuts and jet selection 



For all the results presented here, we consider only a positively charged W and choose 
the leptonic decays, 

W + -> ve + , Z/ 7 * -> e"e + . (8) 

In this paper we shall present results for the Tevatron collider only and we pick a simple set 
of cuts accordingly. All leptons satisfy, 

^ pton > 20 GeV, \y lcpton \ < 1, (9) 

and for the W case there is also a cut on the missing transverse momentum, p™ lss > 20 GeV. 
Our final requirement is that the dilepton mass be greater than 15 GeV. Although this has 
no effect in the W case, it prevents the production of soft e~e + pairs which would otherwise 
be copiously produced by the virtual photon in the Z/'y* process. 



Jets are found using the Run II clustering algorithm [28] with a pseudo-cone of size 
R = 0.7, and are also subject to, 

# > 15 GeV, \y ict \ < 2. (10) 

For the new results on W, Z + 2 jet production, in this paper we will mostly consider events 
where exactly 2 jets are found by the algorithm, i.e. exclusive 2 jet production. The inclusive 
production of jets - which would include events with 3 jets at next-to-leading order - is a 
further option in MCFM that will only be touched on briefly here. 



C. Scale dependence 

The principle motivation for performing a next-to-leading order calculation is to reduce 
the uncertainties in leading order predictions. In particular, any perturbative prediction 
contains an unphysical dependence on renormalization and factorization scales (often cho- 
sen to be equal, as we shall do here). The magnitude of cross-sections and the shape of 
differential distributions can vary greatly between two different choices of scale, which is 
often interpreted as an inherent "theoretical uncertainty" which is then ascribed to the pre- 
dictions. Another strategy is to argue for a particular choice of scale, based on the physics 
of the process under consideration. 
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FIG. 4: The scale dependence of the W and Z+l jet predictions, with the factorization and 
renormalization scales equal and given by /i. The differential distributions dcr/dpx are integrated 
over 15 < px < 200 GeV, with the basic cuts as described in section 1IIB[ 
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FIG. 5: The scale dependence of the W and Z+2 jet differential cross sections, da/dMjj integrated 
over 20 < Mjj < 200 GeV. For more details, see section [II Cj . 
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A next-to-leading order calculation is an invaluable tool for investigating the issue of scale 
dependence. The logarithms that are responsible for the large variations under changes of 
scale at leading order are exactly canceled through to next-to- leading order. As a result, 
one expects that next-to-leading order predictions are more stable under such variations. 
In addition, the next-to-leading order result may provide further evidence to support a 
particular scale choice that may have been deemed appropriate at leading order. 

As an example of expected results, in Figure [§ we show the scale dependence of the 
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exclusive W and Z + 1 jet differential cross-sections da/dpr, integrated over the range 
15 < px < 200 GeV. The next-to-leading order predictions have been known for some 



time [TIJ, but here are calculated within our program, MCFM. For both processes, the 
leading order prediction rises sharply as the scale is decreased, while the corrections produce 
a far flatter curve that exhibits a much less pronounced dependence on the scale choice. 

The corresponding new results for the 2-jet processes are shown in Figure As in the 
1-jet case, the renormalization and factorization scales are set equal and the basic cuts 
mentioned in the previous section are applied. In addition, we now use the dijet mass 
differential distribution, da/dMjj integrated over 20 < Mjj < 200 GeV. As anticipated, 
both processes show a considerable reduction in scale dependence. For example, the ratio 
of the leading order prediction for the W process using a hard scale /J = 2m w to the result 
for a far softer scale /x = mw/2 is, 

cr L o(W + 2 jets, fi = m w /2) = 

( x LO (^ + 2jets,/i = 2m vy ) ' ' 1 ' 

while the same ratio at next-to-leading order is only, 

&nlo(W + 2 jets, fj. = mw/2) _ ^ ^ ^ 
ctnlo(W + 2 jets, \x = 2m w ) 

D. px distributions 

Once again, we repeat some W, Z + 1 jet results, in order to highlight both the similarities 
and the differences with the corresponding 2-jet distributions. 

In Figure |^ we show the jet px distribution for both of the 1-jet cases, using a relatively 
hard choice of scale, \x = 80 GeV. We first note that since we are considering the exclusive 
jet cross-section, the rise of the distributions at low pr is limited only by the jet cut, px > 
15 GeV. At next-to-leading order the distributions change significantly to become much 
softer. At high p? a single jet is much more likely to radiate a soft parton (that passes the 
fixed pt cut and is counted as an extra jet), thus removing it from the sample [ |10| . 

The situation for the 2-jet processes is shown in Figure [7|, where we plot the pt distribution 
of the hardest jet, again using /i = 80 GeV. In contrast with the previous figure, the 
distributions turn over at small px- If the highest px jet has a pt < 20 GeV, there is little 
phase space for the emission of a second softer jet with p T > p™ m = 15 GeV. We also see that 
including the radiative corrections softens the distribution considerably, for the same reason 
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FIG. 6: The jet px distribution of W and Z+l jet events, evaluated with the hard scale choice, 
fi = 80 GeV. 
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FIG. 7: The hardest jet pt distribution in W and Z+2 jet events, using the scale fj, = 80 GeV. 

as before. For the W case, we also show the inclusive distribution, i.e. the cross section for 
the production of two or more jets. This 'fills in' the high-i^r tail of the distribution. 



IV. HEAVY FLAVOR CONTENT OF JETS 



A. Hadron + Hadron -> W/Z + b + b 



We would like to estimate the fraction of W + 2 jet events that contain two heavy 
quark jets. We will limit our discussion to b quarks, because they can be tagged with high 
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FIG. 8: Lowest order diagrams for the process parton + parton — » Wbb. 
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FIG. 9: Sample diagrams for the process parton + parton — > Zbb. 



efficiency. In order to do so, we recall the next-to-leading order results for the production 



of a bb in association with a W, reported in [[L5f| . As a reminder to the reader, we work 
in the approximation in which the b quarks are taken to be massless and we have ignored 
contributions from processes in which there are two b quarks already present in the initial 
state. The basic lowest order diagrams for this process are shown in Fig. |8|. Wbb processes 
accompanied by up to 4 jets have been considered at tree level in ref. |29| . 



For the related study including a Z instead of a W, we use the results presented in Jl6| 



for the production of a bb pair in association with a Z. The same approximations apply as 
discussed above for the W case. The notable difference now is that there are more lowest 
order diagrams, as shown in Fig. 0, including an initial state composed only of gluons. As 



discussed in [|16[], these latter diagrams with initial gluons are believed to be responsible for 
the sizeable corrections to the basic process at large m h i. 

An immediate concern is that neglecting the 6-quark mass may be unjustified J3UJ. At low 
values of m b i, quark mass effects may be important. To address these concerns, in Figure 
we compare the lowest order m h i distribution calculated using the full mass dependence 
with the result obtained by setting nib = 0. There are two effects of introducing a mass for 
the 6-quark. Firstly, the phase-space becomes smaller, leading to a reduction of the cross- 
section. On the other hand, the matrix elements receive extra contributions proportional 
to powers of ml which may increase the result. As shown in Figure [II], the matrix element 
effects dominate around the peak of the distribution where they are as large as 5%. Closer 
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FIG. 10: Illustration of mass effects in lowest order. Mjj is the mass of two tagged 6-jets. 

to threshold the phase space effects are dominant. At large m b i the quark mass effects are 
quite small as expected. 



B. Results 

In order to compare the results for 6-quark jets with those for the whole 2-jet sample, we 
will show the differential cross-section as a function of the dijet mass. We use two choices of 
scale in these analyses, a hard scale \i « M w , M z = 80 GeV and a softer scale fi — 40 GeV. 
These are the scales used for the plots shown in Figures |TT] and |T^|, where both the leading 
order and the radiative corrections are shown for comparison. 

As with the pr distributions of the previous section, the shapes of the distributions change 
when the QCD corrections are included. For both the bb and general 2 jet distributions, the 
hard scale causes the dijet cross-section to increase at next-to- leading order for small values 
of Mjj. The radiative corrections using the soft scale cause a considerable depletion in the 
cross-section at high Mjj. However, the shapes of the bb and 2 jet distributions appear very 
similar when compared at the same order of perturbation theory and using the same scale. 

In Fig. ^ we show the cross section for events that contain 2 fe-tags divided by the 
cross section for all two jet events, as a function of the dijet mass. As can be seen, for 
the W this ratio does not depend very strongly on either the choice of scale or the order 
in perturbation theory. The percentage falls at low values of Mjj until Mjj ~ 60 GeV, 
where it becomes approximately constant at 0.8%. For the Z, the proportion is fairly 
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FIG. 11: Comparison of the expected distributions in the dijet mass for untagged W+2 jet events 
and W+2 6-tagged jet events. A hard scale, \i = 80 GeV, is shown in the left-hand plot and the 
softer scale, \i = 40 GeV, on the right. 
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FIG. 12: Comparison of the expected distributions in the dijet mass for untagged Z+2 jet events 
and Z+2 6-tagged jet events. A hard scale, \x = 80 GeV, is shown in the left-hand plot and the 
softer scale, \x = 40 GeV, on the right. 



constant at approximately 2% for all the curves except for the case of next-to-leading order 
at fi = 40 GeV. In this case the percentage rises at high Mjj. The origin of this effect may 
be associated with the extra diagrams present in the Z case and requires further study. 
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FIG. 13: Ratio of W/Z+2 b-tagged jets to W/Z + 2 jet events in LO and NLO at both fj, = 40 GeV 
and (i = 80 GeV. 

V. CONCLUSIONS 

We have presented the first results for the implementation of W/Z + 2 jet production 
at next-to- leading order in a general purpose Monte Carlo. An analysis based on exclusive 
jet production for Run II of the Tevatron shows that the usual benefits of next-to-leading 
order are realized, among them being a reduced scale dependence and hence an improved 
normalization for distributions. We also find changes in the shapes of distributions similar 
to those found in the 1-jet case. These modifications are reduced if we consider the inclusive 
cross-section. 

We performed an analysis of the heavy flavor content of jets produced in association with 
a vector boson. For production in association with a W, the ratio of 6-tagged to untagged jets 
changes very little upon the inclusion of radiative corrections and appears to be predicted 
very well by perturbation theory. 
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